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Abstract. We propose a new scheme for supplying voltages to the electrodes of 
microfabricated ion traps, enabling access to a regime in which changes to the trapping 
potential are made on timescales much shorter than the period of the secular oscillation 
frequencies of the trapped ions. This opens up possibilities for speeding up the 
transport of ions in segmented ion traps and also provides access to control of multiple 
ions in a string faster than the Coulomb interaction between them. We perform a 
theoretical study of ion transport using these methods in a surface-electrode trap, 
characterizing the precision required for a number of important control parameters. We 
also consider the possibilities and limitations for generating motional state squeezing 
using these techniques, which could be used as a basis for investigations of Gaussian- 
state entanglement. 

PACS numbers: pacs 



Submitted to: New J. Phys. 



Quantum control of the motional states of trapped ions through fast switching 2 
1. Introduction 

Engineered states of the motion of atomic ions in rf traps [U EJ, [3] have provided a number 
of pioneering demonstrations of quantum-state control and decoherence. Examples 
include the creation of Fock-state superpositions, squeezed states, and entangled states 
between the internal and motional degrees of freedom [U El |6j El [9] . The motional 
degrees of freedom are also of primary importance in multi-qubit quantum logic gates, 
where the internal states can be entangled by transient state-dependent excitation of 
the motion [TUl fTT j . More recently, motional degrees of freedom have been used as a 
conduit to transmit quantum states between different trap regions of a multi-zone trap 
array, whether by moving the ions themselves [12] , or by use of the long-range Coulomb 
interaction [13l E] . 

The motion of trapped ions can be controlled by the application of optical or 
microwave fields, or by changing the electric fields which generate the trap itself. 
While optical and microwave fields can be used to create forces which are dependent 
on the ions' internal electronic state (and can therefore be used to produce spin- 
motion entanglement), this control can generally be approximated by a perturbative 
Hamiltonian with respect to that describing the trapping potential. By controlling the 
trap electric fields it is possible to produce large changes in the Hamiltonian describing 
the motion, and this control is independent of the ions' internal state. 

Manipulating the motional states of trapped ions on fast timescales is desirable for a 
number of reasons. The primary application of current interest is quantum information 
processing, where transport of information is necessary for scaling to large numbers 
of ions, and has been predicted to be one of the most time-consuming operations in 
a large-scale processor p~5l [IS]. One promising approach to scaling is to transport 
the ions themselves through an array of microtraps, with the transport controlled by 
time- varying potentials applied to a number of trap electrodes p~3 HE]. These have 
traditionally been generated by voltage supplies placed outside the vacuum system 
and connected to the trap via vacuum feedthroughs. The feasibility of this approach 
has been demonstrated in a number of experiments, including linear transport (and 
separation) of multiple ions [19l [20l [21], as well as 2-dimensional transport through T, 
X and Y junctions [221 1231 EU [25]. The important figures of merit in these experiments 
were the operation time and the amount of motional excitation which persists after the 
transport has taken place. The latter is important because errors in multi-qubit quantum 
logic gates increase for higher motional excitation. Though this can be mitigated by 
sympathetic cooling [26j [271 EH] , it takes time and thus reduces the computing speed of 
the processor [29] . As a result of these considerations, initial experiments on transport 
operated in the adiabatic regime, maintaining the ion close to the ground state of 
motion of the co-moving time-dependent potential well throughout the transport. The 
adiabatic constraint limits transport to timescales slow compared to the motional period 
of oscillation of the ion in the trap. 

The use of low-noise, high-speed Digital to Analog Converters (DACs) has recently 
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enabled diabatic transport of ions over distances of ~ 300 [im in 5 to 16 oscillation cycles 
[201 12T] , with the ion returned to the ground state of the potential well at the end of the 
transport despite being transiently excited during transport. With these techniques, 
the limit to the control rate is set by the finite capacitances present in the lines going 
to the electrodes. 

In this paper, we propose time-dependent control of the trapping potential using 
a new method which involves placing electronic switches inside the vacuum system. 
The control for the switches is digital, allowing trap electrodes to be switched between 
two potentials supplied by standard analog supplies. Since the electronic switches we 
will consider can change their output by up to ^ 10 V on nanosecond timescales, this 
method provides the possibility of making changes to the trapping potentials 100 times 
faster than the period of oscillation of the secular motion, which is typically between 
200 ns and 1 |lls [T7j . This would allow for ion transport in times shorter than a single 
cycle of oscillation in the trap (section [4]). 

Fast control of the trapping potential could also be used for other applications. In 
strings consisting of multiple ions, an additional important timescale is the speed of 
sound in the chain, which is characterized by the frequency splitting of normal modes. 
Controlling the potential on timescales fast compared to the speed of sound may allow 
for the realization of schemes which have been proposed for entanglement generation 
and the investigation of continuous- variable quantum information processing [30l f3Tj . 

This paper proceeds as follows. In section [2] of the paper we introduce our fast 
control scheme. In section [3] we present a generic ion-trap setup which will be used as 
reference for the quantitative studies throughout the text. In section [4] we review the 
application of ultra-fast voltage switching to macroscopic ion transport in less than an 
axial oscillation cycle, while in section [5] we consider a different application to continuous 
quantum- variable manipulation. A number of further applications are outlined in section 
[HJ Finally, a summary of the main results is given in section [7j 

2. Electronic switches for ultra- fast motional control 

The goal of our work is to achieve diabatic ultra-fast (for our purposes, this is defined 
as being faster than 100 MHz) control of the trap potential of one or more trapped 
ions. Additionally, we do not want to relax the requirements on voltage noise on the 
electrodes, which can lead to ion heating, setting a limit to the coherence of the quantum 
states of motion and shortening the lifetime of the ion in the trap. The rate of excitation 
of a single ion from the ground state to the first excited state (which characterizes the 
rate of heating in a trap) is given by [32] 

^2 



where q is the ion charge, m its mass and uu is the secular frequency in angular units. 
Se(w) = 2 /f^ dr (E(t)E(t + r)) e luJT is the spectral density of electric field fluctuations 
in the trap, where t is time and (• • •) time-averages the argument. In order to prevent 
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Figure 1. Electronic schemes for setting the voltages at the control electrodes of a 
trap. The standard configuration is shown in a) and consists of a voltage supply and 
a low-frequency filter (rm, a and Cfn, a can be chosen such that the corner frequency is 
f c ~ 10 Hz). The configuration in b) allows for waveforms, typically at frequencies 
close to but lower than the secular frequency. The idea is the same as in the standard 
configuration, only the voltage supply is replaced by a digit al-to- analog converter and 
the corner frequency of the low-pass filter is set to a higher value (f c ~ 100 kHz). The 
ultra-fast switching scheme described in the text is shown in c). The values of r^ c 
and Cfii ?c can be chosen such that f c ~ 10 Hz. In all cases, the capacitors have two 
functions: they form part of the low-pass filters for preventing high-frequency noise 
from getting into the trap; and they serve as a low-impedance path for the rf drive 
to ground. c tra p is the capacitance between the rf and the control electrodes, and 
Ccabie represents the parasitic capacitance of the typically long cables and feedthrough 
assemblies from the external electronics to the trap chamber. The small insets next to 
the supplies indicate the generic time-dependent voltages on the respective lines. 



relevant noise spectral components (hundreds of kHz to tens of MHz) from reaching a 
trapped ion, it is common to place low-pass RC filters close to the ion trap (see figure [I]). 
These filters attenuate any contribution to Se{u>) which comes from outside the trap. 
However, they also impede the ability to perform fast control of the trapping potentials 
from outside the vacuum system, since they suppress high frequency components in the 
analog signals and therefore distort them. 

One solution could be to replace the low-pass filters by band-stop filters which allow 
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Figure 2. Switching behavior of a CMOS SPDT switch (74HC4066M from Texas 
Instruments) at 4 K. The 10%-90% rise time of the switched output is below 4.0 ns. 
The upper part of the curve is the expected response due to the load of the detection 
electronics (the coaxial cables and the oscilloscope introduce a capacitive load of 
~ 50 pF), which would not be present in the wiring scheme of a trap electrode. 



for both slow and very fast signals. However, parasitic capacitances of the cables to the 
electrodes (including vacuum feedthroughs) also filter the higher-frequency components 
and inhibit ultra-fast control. Another possibility is to adjust the waveforms to the 
filters, such that the voltage control sources deliver high amplitudes of the fast transients 
and a fraction of them reaches the electrode [33]. In this case, the limit is set by the 
compromise between the filter strength and the maximum amplitude, and the slew rate 
of the voltage source. 

As an alternative to these approaches, we propose to generate the time-dependent 
control close to the electrodes themselves. The primary ingredient is a single-pole 
double-throw (SPDT) semiconductor switch, the state of which is controlled by a digital 
line. The SPDT switch provides a connection from one of two input lines to the output 
dependent on the logical state of the digital input. The chosen input switches abruptly 
when the digital input reaches a threshold value, thus the precise details of the digital 
signal arriving at the switch are not important as long as the time at which it reaches 
the switching value does not fluctuate. The speed at which the voltage switches is 
determined exclusively by the characteristics of the switch itself and the capacitive load 
at the output. 

Figure [j] c) shows the electronic scheme for wiring up a control electrode to analog 
inputs in our proposal. The switch is placed between the low-pass filter and the 
electrode, otherwise the filter would also attenuate the fast components of the switched 
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signal. 

An illustration of the analog output of a candidate CMOS SPDT switch 
(74HC4066M chip from Texas Instruments) operating at 4 K is shown in figure [2j taken 
from tests of a range of technologies which will be published elsewhere p3]. The digital 
control signal is not shown, but it is much slower (~ 50 ns). Despite a slower the digital 
signal, the output switches in ^ 4 ns. Also visible in the signal is a slow component 
which changes over a time of 15 ns. This is due to our measurement setup, including 
the capacitance of a 1 m co-axial cable which was used to connect the output of the 
switch to the high-impedance oscilloscope and other electronics, introducing an extra 
capacitive load (~ 50 pF). 

As a result of our proposed wiring scheme, the route from the electrode to rf ground 
now includes the internal impedance of the switch on the on state. It is important that 
this impedance is low to pull the electrode to rf ground. The on-resistance of the 
74HC4066M is specified to be ^ 15 Q at 300 K [35], and is probably smaller at 4 K, 
for which we found the switching time to be 3 times faster than at room temperature. 
As a reference, the typical reactance of the capacitors in the RC-filters ranges from 
0.5 — 10 for trapping-drive frequencies of 20 — 200 MHz and capacitors close to 1 nF 
(above ^ 1 nF parasitic inductive effects are relevant). These values have been used in 
previous experiments [36j. 

Another relevant aspect is to design the electronics such that the coupling between 
the digital and analog sides of electronics is minimized and noise on the digital lines 
does not reach the electrode (cross talk). 

3. Experimental implementation 

As an illustration of the considerations which must be taken into account when working 
in this new regime, we next consider the application of our scheme to two tasks in 
quantum control: diabatic transport of ions and the generation of squeezed states of 
the motional degrees of freedom of ions. In order to provide an experimental setting, 
we introduce the fixed geometry of a surface-electrode ion-trap [371 EH [39] which we 
have recently fabricated. This trap design provides a setting which gives an idea of the 
size of relevant experimental parameters. However, the methods and routines which we 
give are general and can in principle be carried out in any other ion-trap setup, so long 
as there are a sufficient number of independently controllable electrodes. The central 
region of the trap is shown in figure [3| The rf-electrodes and the control electrode e m 
run parallel to the trap axis. By applying an rf voltage with an amplitude of 100 V 
at a frequency of 100 MHz, this structure yields a three-dimensional confining pseudo- 
potential for positively charged particles at a trap-electrode distance of ~ 60 \im. The 
array of control electrodes allows for flexible creation of potential wells along the axis 
of the trap. 
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Figure 3. Surface-electrode trap suitable for ultra- fast quantum-control experiments. 
The rf-electrodes are labeled rf and the control electrodes e. The main trapping zone 
is the area sandwiched between the five pairs of narrower (50 |Ltm) control electrodes, 
denoted as e2-8- 



4. Fast ion transport 

The basic principle of fast transport using sudden switching of potentials is illustrated 
in figure [4] for an ion considered a 1-dimensional harmonic oscillator [40i I4T] . The ion 
starts in the ground state of an initial potential well situated at z = — z given by 

VuAt(z) = ^mu; 2 (z + Zo) 2 . (2) 

At time t — the potential is suddenly displaced to the transport well, which has the 
same curvature but is centered at z = 0. At this point the ion is in a coherent state of 
the transport well, with coherent state parameter a = — z /a , where a = ^Ti/2muj is 
the r.m.s. extent of the ground-state wavefunction. Under free evolution, this coherent 
state will gain and lose momentum, returning to rest periodically at times t = pn/uj 
where p is an integer. For p — 1, the wavepacket is positioned at z = Zq. If at this time 
the potential is suddenly displaced again to a final potential well which is centered at 
z = Zq and has the same curvature as the transport well, the ion will end up in the 
ground state of the final potential, having been transported over a distance of 2z . 

Such a "throw-catch" transport routine can be carried out in the trap shown in 
figure [3j One possibility is to start with a potential well centered above the pair of 
electrodes e 6 and end above the pair e 4 . The full transport distance would then be 
100 |xm. At a typical frequency of 1 MHz, this gives a mean velocity of 200 m/s and 
a - 4500 for a 40 Ca + ion. 
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Figure 4. Sketch of the fast "throw-catch" transport routine, (a) Starting point: the 
ion is in the motional ground state of a potential well displaced from the trap center 
by an offset voltage on an external electrode, (b) Ion transport: the offset voltage is 
"instantaneously" switched off so that the potential is centered in the trap and the ion 
starts its coherent oscillation, (c) Ion catch: the offset is switched onto a symmetric 
external electrode and the ion is caught at the instant when it has no kinetic energy, 
exactly half an oscillation after beginning the transport. 

There are a number of experimental challenges behind this conceptually simple 
transport routine, which we outline below. 



4-1- Micromotion 

Micromotion is the driven motion of trapped ions in rf traps which occurs when the rf 
component of the electric field is non-zero. 

Radial micromotion arises when the dc minimum does not lie on the rf minimum, 
which is given solely by the geometry of the trap. There are two causes of 
radial micromotion: stray fields and trapping dc voltages. The former are a priori 
unpredictable and should be canceled experimentally. In the transport routine, the effect 
of changing the potential of any single electrode is not only to displace the ion along the 
z-axis, but also in the radial directions. Thus, it will be necessary to switch multiple 
electrodes simultaneously both in steps (a) and (c) in figure [ij The set of voltages 
required in order to cancel radial micromotion can be determined from simulations 
[421 [431 144] . As an example, we give in table [j] the results of such simulations for the 
trap described in section [3) 

Axial micromotion arises when the axial component of the rf-drive does not vanish, 
usually due to misalignment, manufacturing imperfections or finite electrode structures. 
There are two relevant components to axial micromotion: one arising from a residual 
electric field of homogeneous amplitude along the axis, and the second due to the 
curvature of the pseudo-potential along the trap axis. For an ion of charge Q and mass 
m, the combined effect leads to a time-dependent motion which can be approximated 

by 145! 

x(t) ft C cos(cj£) + C 2 cos[(fi + u)t] + C-2 cos[(f2 - u)t] + D 2 cos(ft£),(3) 
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Table 1. Trap voltages required to implement the fast transport routine in the trap 
shown in figure [5] with a single 40 Ca + ion. The transport takes place from e§ to 
(100 |Ltm) in half an axial oscillation (500 ns). The largest voltage switch is < 6 V. 



The energy E of the ion and the relevant potential parameters (defined in section 4.2) 
are also given. Some electrodes require three different voltages, for which two SPDT 
switches can be used. Note that the numbers are particular to a specific trap geometry 
and conditions, but may be used as a reference for different scenarios. The labeling of 
the electrodes is given in figure [3] 
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where Q and uj are the drive and secular angular frequencies, respectively, Co is the 
(classical) amplitude of the secular oscillation (100 |Lim in the transport example above), 
C±2 are due to the curvature of the pseudo-potential and D 2 is due to the residual 
oscillating electric field. The size of C± 2 is given by 

c *"- s,-frW <4) 

with 



f3 ny/a z + q*/2, 



SQ 
4Q 



mQ? 2 ' 



Here, c 2 and d 2 are the second order coefficients in the dc and rf contributions to the 
axial potential, respectively. In our trap, c 2 — 8 • 10 -6 V/m 2 and d 2 — 4V/m 2 , so 
C± 2 ~ 2.5 fm. D 2 could be larger, but its effect will time-average to zero during ion 
transport, since the drive frequency is required to be much higher than the axial secular 
frequency for stable trapping. 

In practice, a careful experimental characterization of the potentials created by the 
electrodes will be required in situ before the proposed methods can be implemented, 
because any real electrode structure will be finite, have gaps and be subject to fabrication 
imperfections. 

4.2. Effect of trap anharmonicity 

For macroscopic transport distances (100 |xm in the example considered), the 
anharmonicities seen by the ion during its coherent oscillation might not be negligible. 
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Table 2. Voltage combination which cancels even terms in the Taylor expansion of 
the axial potential up to 10 th order. Note that this is a particular solution for which 
the voltages at e^i are constrained to be equal to those at e- hr for all electrode pairs 
except for e 4 _6- All voltages are given in volts. The labeling of the electrodes is given 
in figure [3j 
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The trapping potential can be Taylor-expanded as 

V(z) » 2n 2 mf*z 2 (l + + sgn(L 4 )|^ , (6) 

where f z is the oscillation frequency of the ion (determined from the curvature of the 
potential), and L 3 and L 4 are, respectively, the length scales at which the effects of the 
cubic and quart ic terms of the potential are of the same size as the quadratic term. 
The principal anharmonicity relevant in the dynamics described above is a quartic term 
which causes a variation of the curvature of the potential between z = ±z and z = 
(the odd term in the expansion in equation [6] is negligible due to symmetry, see table 
[I]). The effect of L 4 on the frequency can be evaluated by taking the second spatial 
derivative of the potential in equation [6| Simple calculations lead to 

^)=^(0)(l + 6^^). (7) 

The consequences of equation [7] were simulated numerically using a Suzuki- Trotter 
expansion [47J of the full Hamiltonian (split-operator method). Due to the size of 
the coherent state under consideration, the simulation required a time resolution of 
360 ps (1400 time steps) and a position resolution of ^ 8 pm to describe the wavepacket 
(r.m.s. extent of ~ 10 nm). 

Figure [5] shows the effect of the quartic anharmonicity on the probability of overlap 
of the final state with the lowest Fock states in the final well. All three potential 
wells were taken to have the same curvature at their respective centers, leading to an 
axial frequency of 1 MHz. For negative values of L 4 , the curvature of the transport 
potential becomes negative at \zi\ > |L 4 1/^/6 and the potential becomes anti-confining 
at \z2 1 > |L 4 |/\/2- In general, the final state is no longer a minimum uncertainty state, 
since the variance in momentum increases. Using the parameters listed in table [TJ the 
probability of overlap of the state in the catch potential with its ground state is 90%, 
and with the n = 2 state 10% at t = 500 ns. 

One possibility to minimize this effect is to design the transport well such that 
the quartic term in the potential is reduced, either using geometry or suitably chosen 
electrode potentials From simulations we find a combination of voltages for which 
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Figure 5. Probability of overlap of the final state with the Fock states n = 0,2,4,6 
of the final potential as a function of the strength of the quartic anharmonic term 
during transport. The data points have been simulated for a 40 Ca + ion at an axial 
frequency of 1 MHz which undergoes a transport of 100 |Ltm (zq = 50 |Ltm). For 
1/L 4 < —0.014 |Ltm _1 , the potential becomes anti-confining at the edges of the 
transport region. 

all contributions up to 10 th order are canceled, leading to frequency shifts below 3% 
between z = ±z and z — 0. However, it would require applying potentials spanning 
~ 80 V, see table [2j This is currently not possible due to technical limitations of 
the switches. In addition, it requires further electrodes to compensate for the radial 
distortion of the pseudo-potential. 

We also investigated having a larger curvature in the initial and final wells than in 
the transport well. In that case, the ground state of the initial well projects onto 
a squeezed (rather than coherent) state of the transport well. The probability of 
overlap with the ground state of the final well is shown in figure [6] as a function of 
the frequency of the initial and final potentials, given a fixed transport potential. The 
highest frequency simulated was 3.0 MHz due to computational limitations. Until that 
value, the probability of overlap increases because the effect of the anharmonicity of 
the transport potential is smaller across the squeezed wave function. It might come as 
surprising that the curve seems to asymptotically approach a probability of overlap of 
100%. Indeed, we expect that it will reach maximum. The reason is that, close to the 
trap center, the wavefunction of the squeezed state broadens. Above a certain frequency, 
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Figure 6. Probability of overlap with the ground state of the final potential as a 
function of the frequency of the initial and final potentials. The transport potential is 
fixed so that u;(0) = 2tt • (1 MHz) and L4 = —120 |im. The total transport distance 
is 100 |Lim (zq = 50 |im). The highest frequency simulated was 3.0 MHz due to 
computational limitations. 



the effect of the anharmonicity on the broadened wavefunction will be large enough to 
reduce the probability of overlap with the final potential well. 

4.3. Switch timing 

Since the timing of the switch to the final potential well is a critical part of the transport 
protocol, it is worth examining the tolerance of the scheme to timing imprecision. We 
can make an estimate by considering that the ion in the transport potential is in a 
coherent state characterized by the eigenvalue 

ao = -- (8) 
a 

of the harmonic oscillator annihilation operator, where Zq is the initial displacement 
from the midpoint and a = \J%j {2muj) as before. If the catch potential is switched on 
at if = tt/u) + St (rather than at t = tt/uj), the ion was in the coherent state 

\a { ) = \-a e-^ St ) (9) 

of the transport potential. The overlap with the ground state of the final well is then 
given by 

<0| D\-a ) \a t ) = (0 |-a (-1 + e~ iujSt ) ) ~ exp (a uj5t) 2 } . (10) 
where we have approximated St <C 1/uo. 
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Figure 7. Position of the potential minimum as a function of time during transport. 
The transitions from to r and from T to T + r' will be given by the experimental 
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The probability of overlap is Po = (0| D*(—ao) {a?) , so the constraint on St in 
order to achieve an overlap of Po is 

St = ^^. (11) 

For a 40 Ca + ion at u = 2n • (1 MHz) and Zq = 50 |xm, the r.m.s. extent of the initial 
ground-state wavefunction is a — 11 nm, the initial coherent state is \a \ ~ 4450 and the 
timing resolution required is St < 12 ps for Po > 0.9. This poses a demanding constraint 
on the trapping potential, which should be stable to ~ 10 ps/1 |Lis = 10 -5 throughout a 
measurement to avoid recalibration of the digital delays within an experiment. Timing 
resolutions of the order required are achievable with commercially available digital delay 
generators j48j [49] . However it should be noticed that delays due to impedances in the 
connections to the trap and the trap itself must also be stable at this level so that they 
can be calibrated out. 



4-4- Finite switching time 

Until now, we have assumed that the voltage switching is instantaneous. In our proposed 
transport scheme, the ion is transported by two subsequent displacements of a constant- 
curvature harmonic potential, which we can describe by \muj 2 {z — s{t)) 2 , where s(t) is 
the position of the potential minimum. Finite switching times of electrode potentials 
will mean that the potential minimum takes a certain time r to move from the initial 
position s(0) = — z to the transport well centered at s(r) — 0. The same applies to the 
catch procedure, which will start at s(T) = and will end at s(T + r f ) = z (see figure 
[7]). Note that r and r f are related to the switching time of the electrode potentials, but 
may slightly differ from it if more than one electrode is switched. 
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Figure 8. Black solid line (left axis): minimum size of the final coherent state \a\ as a 
function of r' . Red dashed line (right axis): probability of overlap of the final coherent 
state with the ground state of the final potential. Data calculated for a 40 Ca + ion 
at uj = 2ir • (1 MHz), zq = 50 |Lim and r = 5 ns. The displacements of the potential 
minima are assumed to be linear. 



In [50] . it is shown that if the ion is initially in the ground state of motion on the 
starting potential s(0), it will evolve into a coherent state with amplitude 

rt 



a(t) = - 



muj 

~2h ( 



-iujt 



dt's(t')e 



iujt 



(12) 



with respect to the frame centered at s{t). For our case, this means that the ion will 
finish in a coherent state with amplitude 



a(T + r') 



~2h ( 



-IL0(T+T r ) 



+ 



r' 

/ dt's(t' + T)e 
Jo 



dt's(t')e iLOt 

ioj(t'+T) 



(13) 



with respect to the frame centered at s(T + r f ). 

From equation [13] it can be seen that the effects of finite switching times can be 
canceled out, i.e. it is possible to catch the ion in the ground state of a final potential, at 
the expense of having T and s(T + r f ) as variable parameters. However, it is in general 
not possible to make the transport symmetric, i.e. fix s(T + r f ) = — s(0), and at the 
same time catch the ion in the ground state by controlling only T. 

In order to estimate the size of the effects of finite switching times, we performed 
analytical calculations (see Appendix A) for the case of linear and sinusoidal transitions. 
For realistic experimental parameters their effect is very similar. As an illustration, the 
size of the final coherent state for the linear case (equation A.l) is shown as a function 
of r f in figure for a 40 Ca + ion at uj = 2tt • (1 MHz), z = 50 |am and r = 5 ns. In such 
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Figure 9. a) Connections required to switch multiple trapping regions, in an ideal case. 
One row of electrodes and connections to the bus are shown. The ions are separated by 
three electrodes in this scheme. Electrodes A are switched between voltages H and LL, 
and electrodes B between and L, with H>0>L>LL. Control lines Ca and Cb switch 
the voltages applied to the pairs of electrodes A and B, respectively, b) Switching 
steps for the transport operations. The ions' positions (centered above an electrode) 
are indicated by a shaded cell. It allows an ion to be moved two electrodes to the right 
with a time period tt/uj. Ati account for the time delays during which the ions are 
in the ground state of the potential wells centered above the B electrodes. AU = 
would be desirable for fast transport, while finite values would allow for operations to 
be carried out on the ions. 



a setup, to catch the ion with a probability of overlap above 90% with the ground-state 
wavefunction, it suffices to have r f < 1.5r. 

For an arbitrarily shaped transition, if s(T + r') = Zq, T = tt/uj and r' = r, then 



equation 13 becomes 



a(T + t) = ]f^e- ZUJT £ df [s(t') - s(t' + T)] e^' . (14) 

Hence, if the transitions are equivalent, i.e. s(t') = s(t' + T), we can recover the ion in 
the motional ground state after transport, making the process first-order insensitive to 
finite switching times. 

4-5. Extended transport throughout an array 



In a large-scale trapped-ion quantum information processor based on the architecture 
described by Kielpinski et al [18], it is necessary to move ions between many remote 
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trapping regions. Therefore it is worth considering how the scheme outlined above could 
be extended and implemented in a repeated fashion. Here, experimental simplicity 
means that it is desirable to minimize the number of connections which are required. 
In figure [9] we give a simplified outline of a switch fabric which could be used to wire 
up electrodes in order to implement transport between a number of adjacent zones 
of a multi-zone trap. Though this does not include all electrodes required to control 
transport over any single zone, it illustrates that increasing the number of zones does not 
imply a direct increase in the number of digital and analog control lines, since repeated 
transport in this scheme is realized by a total of four analog voltages controlled by three 
digital lines. 

In a real experiment, it is likely that fabrication imperfections and stray 
environmental electric fields would require additional electrodes which could be used 
to tune the trap frequency of the ion. In addition, switch timing could be controlled by 
adding phase-shifting elements such as varactors, which have previously been proven to 
work in cryogenic ion trap settings [53] . 

With a chain of ions of equal mass, it should be possible to shunt the whole string 
along by one trap zone by switching the electrodes in the same manner as for a single 
ion. This is due to the fact that, as long as the curvature of the potential seen by 
the ions during the transport process remains the same, no motional modes will be 
affected other than the center-of-mass mode, which is excited exactly as a single ion 
would be [52]. In addition to transport, in a large-scale processor it is likely that the 
ability to separate and deterministically re-order ion strings will be necessary. It is 
conceivable that methods similar to those we have described above could be used for 
these tasks. Separation and recombination will inevitably involve anharmonicity in the 
potential (trapping and Coulomb repulsion), which may require methods such as those 
we describe in section [42] to keep the ions close to the ground state. 



5. Motional-state squeezing 

An alternative application of the fast switching method is in the generation of non- 
classical states. Previous work by Serafini and co-workers describes methods for 
generating continuous- variable entanglement by applying ultra-fast control of the radial 
degrees of freedom of an ion chain [301 [31] . A critical ingredient of the Serafini scheme 
is the generation of squeezing of one of the ions, with the Coulomb interaction between 
ions then giving rise to entanglement. As a first step towards such experiments, we 
consider the lifetime of the states squeezed by applying ultra-fast control of the trapping 
potential, which is limited by the presence of a quartic term in the potential as well as 
heating due to fluctuating electric fields. 

Squeezing refers to a state in which quantum fluctuations in one quadrature 
are suppressed relative to their value in the ground state [54] . The orthogonal 
quadrature has increased fluctuations, preserving the Heisenberg uncertainty relation. 
Mathematically, a squeezed state can be obtained from the vacuum state by applying 
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Figure 10. Illustration of the squeezing mechanism based on voltage switching. The 
ion starts in the motional ground state of the first potential and at t = its curvature 
is suddenly switched. The wave function starts breathing periodically with angular 
frequency 2co>2, accommodating to the new potential. At t = tt/ (2,002) it reaches a 
maximum width and the potential is switched back, leaving the state squeezed. Black 
solid lines represent the wave function and the potential at the beginning of the period 
indicated and red dashed lines at the end of it. 



the operator 

S(0 = e (^) 2 -^ f ) 2 )/ 2 , (15) 

where the parameter £ = re l<t> is the complex squeezing parameter, denoted by a real 
magnitude r and phase <£>. The squeezed state is one example of a Gaussian continuous 
variable state, which can be characterized in terms of the first moments TrJTt^-p], and 
the second moments, embodied in the covariance matrix a 
1 



-Tr[{^-, R k }p] - Tr[R jP ]Tr[R k p}. 



(16) 



Here, p is the Gaussian state and R is a vector of normalized position and momentum 
operators. For one ion, R = (x,p) T . The amount of squeezing of the state is given by 
the smallest eigenvalue of the covariance matrix. For details, see e.g. [57] . 

Squeezing of motional states of trapped ions has already been reported using optical 
fields to create a parametric drive which squeezed the vacuum state [4]. Instead, we 
consider the squeezing procedure proposed in [30]. It involves an abrupt change of the 
trapping frequency from u to \uj at time t — 0, followed by the opposite transformation 
at time t = h/(2\uj) (see figure 10). For a Gaussian initial state with mean zero, the 



time evolution can be fully characterized by the transformation of the covariance matrix 
according to a t = S(t)aS(t) T , where [30] 



S(t) = 



\ 



\ 



cos(oj\t) A sin(cjAt) 
—A -1 sin(cc;A£) cos(cc;Ai) j 
is a symplectic transformation. For t = h/(2\uj) we find 

/ X2 \ 

A (T22 — 021 

Ot = 



(17) 



A (T22 — 021 
-012 01l/A 2 



(18) 



Quantum control of the motional states of trapped ions through fast switching 18 



i.e. the variance of one quadrature is squeezed by a factor A 2 if A < 1, while the other 
increases by the inverse amount. 

The major and minor quadrature normalized variances can be shown to be e 2r h/2 
and e~ 2r Ti/2 ) respectively [58J. Therefore, the maximum r.m.s. extent of the squeezed- 
state wavefunction for a trapped ion is simply Ax = a e r = a /A, where the last equality 



comes from equation 18 



The trap frequency is proportional to the square root of the voltage applied 
to the central trap electrodes (e 5 in figure [3]), given that the rest of electrodes are 
grounded. For the trap described above, switching the trapping voltage from 10 V 
(at! ~ 2tt • (1.2 MHz)) to 1 V (uo 2 ^ 2tt • (400 kHz)), waiting for a time of tt/(2cj 2 ) ^ 
625 ns and switching back to would give A = \/l0 ~ 0.32 vacuum units or —5 dB 
(A = dB corresponds to the non-squeezed case). Furthermore, the wavefunction of 
a squeezed state breathes at 2cji, so the squeezing procedure described above can be 
repeated with a period At = 71-/(2^2) + tt/(2uji) = (1 + 1 / X)tt / (2uji) to enhance the 
results. 

In the following, we consider the two main experimental limitations to such a 
procedure: motional heating of the ion in the trap, and anharmonicities in the trapping 
potential. 

5.1. Motional heating 

An important factor which limits quantum-state engineering of the motional states of 
trapped ions is heating due to fluctuating electric fields at the position of the ion [59, 60j. 
In order to estimate the degree to which heating might limit squeezing using the methods 
discussed in the previous section, we consider the effect of a time-varying electric field 



E(t) on a squeezed state. In Appendix B| (equation B.16) we show that the loss of 



fidelity with the original squeezed state is characterized by a time constant r, which we 
can rewrite as 

_ x cosh(2r) e 2r T ^ 

r = IW— 3— « IW T = — , (19) 

where r ^i is given in equation [T] and the approximation is valid for r 3> 1. 

Extrapolating heating-rate measurements in [61] to the temperature and ion- 
electrode distance for the trap in figure [3] yields heating rates down to r ^i ~ 
10 quanta/s. For A ~ — 10 dB, which is the maximum squeezing obtained with other 



continuous-quantum- variable systems to date [62], equation 19 yields a lifetime of the 
squeezed state of ~ 4 ms, much longer than the creation time for this state in our 
scheme. 

5.2. Effect of anharmonicity 

Anharmonicity in the trapping potential presents another limit to the amount of 
squeezing which could be achieved. Finite quart ic terms in the potential are the most 
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relevant anharmonic contribution for a symmetric trap (see table [T]) and result in an 
effective spatial frequency distribution cj(z) as given in equation [7J 

In the phase-space representation given by the Wigner distribution [63] , a squeezed 
vacuum state takes a two-dimensional Gaussian form centered at the origin [6lj. If we 
let it evolve in time in a harmonic potential, it will rotate around the origin at the 
trap angular frequency uj. However, in an anharmonic potential the wave function will 
distort. To quantify the lifetime of the squeezed state in the anharmonic trap, we can 
calculate the overlap between its wave function and the one it would have become had 
it evolved in a harmonic potential: 



Po(t) = 



2 



(20) 



We have used the same numerical techniques as in section |4.2| to carry out 
simulations for the lifetimes of squeezed states. For a 40 Ca + ion at uj = 2tt - (1 MHz) 
with A = — 10 dB and L 4 = —120 |xm, the lifetime is ^ 3 ms. This is comparable to the 
lifetime expected from the motional heating effects considered in the previous section. 



However, as discussed in section |4.2[ it is possible to improve on the harmonicity of 
the trapping potential by optimizing the trap geometry and the trapping potentials, 
whereas reducing motional heating effects is experimentally very challenging. Therefore, 
we expect that heating effects, not anharmonicities in the potential, will ultimately set 
the limit to the lifetime of the squeezed motional states in the trap. 

6. Outlook 

In the following, we outline a range of further applications which could profit from 
ultra-fast motional control. 

6.1. Entanglement between radial motional modes 

Linear couplings between two harmonic oscillator modes which are prepared in Gaussian 
states squeezed along orthogonal quadratures can give rise to entanglement [56l |65] . In 
the case of a string of trapped ions, the motional oscillations of the individual ions are 
coupled by the Coulomb interaction, for which the lowest order term in the expansion 
of the interaction potential energy is 

A* 6 *- (2i » 

where £i and £ 2 are the excursions of the ion about equilibrium in a given direction, and 
X = 2 if the oscillation is aligned with the vector connecting the equilibrium positions of 
the two ions and \ = 1 otherwise. For two identical ions trapped in the same potential 
well, the distance between the ions is 

^(adb) 1 * (22) 

where m is the mass of a single ion. 
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Entanglement can be characterized when the two subsystems which are entangled 
can be accessed locally on timescales which are fast compared to the interaction between 
these systems. For the case of ions trapped in a single potential well, the interaction 
energy can be characterized by the splitting of the normal modes, which is typically 
on the order of a few hundred kHz, and thus the local access to variables (in this case 
changes in the potential) must happen on timescales short compared to 1 |Lts. Once the 
voltages are set back to the initial configuration, the squeezing will be transferred from 
local quadratures to a non-local combination of the coordinates, which will lead to the 
entanglement. 

6.2. Dynamical studies 

A trapped ion is subject to "kicks" due to elastic collisions with background-gas 
atoms and molecules. In order to understand the mechanisms that lead to Paul-trap 
instabilities, it might be interesting to map the effect of such kicks as a function of the 
collisional energy and the resulting impulse direction of the ion. 

Collisions could be emulated by pulsing on electric fields. The effect of such 
collisions can be created in our trap with a routine that starts by cooling an ion to 
the motional ground state and then adding a sudden offset potential to one or more 
electrodes, in order to displace the potential minimum and kick the ion in a certain 
direction. If the ion is left to evolve in the new potential for a time 7r/(2a;), it will 
acquire the maximum kinetic energy. At that moment, suddenly removing the offset 
potential would leave the ion with a high kinetic energy in the original potential well. 
Another possibility would be to let the ion evolve for a time tt/uj in the second well, and 
then remove the offset to leave the ion with a higher potential energy in the final well. 

At room-temperature, the mean thermal energy of the background-gas particles is 
E t h ~ 25 meV, so the offset voltages required for the routine above are very similar to 
those required for the ion transport routine (voltages below 10 V on electrodes e 2 -8 are 
enough for energy kicks above 20 meV, see table [j]). 

6.3. Interface with solid-state devices 

Fast voltage switching could also be used to help in realizing an interface between a 
trapped ion and a solid-state quantum device, which could be used to map quantum 
information from a solid-state qubit to an ion and vice- versa [66] . In order to enhance the 
coupling, it is desirable to place the ion close to the solid-state device, where anomalous 
heating could be significant. A fast transport routine similar to that described in section 
[4] could be of use to quickly bring the ion close to and away from the device, thus limiting 
heating effects which occur outside the desired interaction time. 
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7. Conclusion 

We have proposed a method for controlling the potentials in a radio-frequency ion trap 
on nanosecond timescales. This would allow for ultra- fast control of the potentials seen 
by the trapped ions, i.e. control at rates much faster than the ions' secular oscillation 
frequencies. The switching time in our proposal is limited by the capacitive load on the 
semiconductor-based voltage switches used. We have experimentally measured this to 
be around 4 ns in a 4 Kelvin test apparatus with a load of ~ 50 pF. 

One possibility opened up by this new approach is to transport an ion within a 
single oscillation cycle. We have shown that this is experimentally challenging but 
should be possible for transporting a single 40 Ca + ion at an axial angular frequency 
uj = 2tt • (1 MHz) over a distance of 100 |xm in 500 ns. 

Another application is the generation of motional squeezed states, where squeezing 
the ground wavefunction of a 40 Ca + ion at uj = 2tt • (1 MHz) by a factor A = —10 dB 
should be within experimental reach. 
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Appendix A. Finite switching time - linear and sinusoidal potential 
displacements 

Here we calculate and then compare the effect of two different transitions between the 
initial and transport potential wells for the case where voltage switching takes a finite 
time (see section 4.4). 

Let us consider first that the transitions from t = to r and from t = T to T + r f 
are linear (as in figure [7]), such that their slopes are Sthrow = W r an d ^catch = ^o/ 7 " 7 ? 
respectively. In this case, equation p~3| becomes 

a lin (T + r') = 



/ lin (T) 



UJT 



^-e-^7 Un (T), 

2TlUJ T 



UJT\ 



UJT' 



2 J 



(A.l) 



where the approximation is for r, r' <C 1/uj. The value of T which minimizes f lm (T) 
has the analytical form 



-ilin 



Illl <l log 

UJ 



2i — UJT 
-2i + ujt' 



(A.2) 



and yields an |a^ n (T + r')| shown on figure fel for the specific case of a °Ca ion at 



UJ 



2n • (1 MHz), zq = 50 |xm and r = 5 ns 
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A sinusoidal transition, 

^throw(^) 



^catch(^) 



z 



2 
1 



1 



/Tit 

cos b 



(A.3) 



resembles more closely the pulse shown in figure [2| while still yielding an analytic 
solution: 



of 



5 (T + r') 



mu itzq^. ecos 



f cos (T) 



2h 2 



IUJT 



/ cos (T), 



2 _ e -iu T + 



IUJT 



7T 2 — (OUT) 2 7T 2 — (bJT f ) 2 

where the approximation is again for r, r' <C 1/ui. In this case, 



1 



Im < log 



(7T 2 -UJ 2 T ,2 )(-4-2lUJT + UJ 2 T 2 ) 



(A.4) 



(A.5) 



u I - v (tt 2 - cj 2 t 2 ) (-4 - 2iur f + uj 2 t' 2 ) J J 
A comparison of the amplitude of the final coherent states for the two different 



examples (equations A.l and A.4) is shown on figure Al, for the same conditions as 
above. The small difference between them seems to indicate that the exact shape of 
s C atch,throw(^) is not greatly important. 



Appendix B. Heating rate of squeezed vacuum states 

Here we will calculate the lifetime of squeezed vacuum states, which can be generated 
in a trapped-ion experiment as discussed in section [5j We will start by considering the 
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transition rate out of an arbitrary state |</>), given by 

|(|(0I^)I0)| 2 ), (B.I) 

where U (t) is the time-evolution operator and the right and left angle brackets average 
over all realizations. 

We consider here the effect of a time-dependent perturbation Hamiltonian H(t) 
which can be written as 

H(t) = -qE(t)z, (B.2) 

where q is the charge of the ion, E(t) the electric field, and z is the coordinate of the 
ion center of mass. We restrict our study to a single dimension, since the normal modes 
of a trapped ion are independent and thus the extension to other dimensions is trivial. 
Heating (or decoherence) is measured as an average over many runs of an experiment, so 
E(t) is a stochastic variable. If it is derived from a stationary process with fluctuations 
around a mean value of zero, we can write averaged quantities 

(E(t)) = (B.3) 
(E(t)E(t + T)) = C(r), (B.4) 

where C is a function to be determined experimentally, but which depends only on r. 
A Hamiltonian of the form given in equation |B.2| gives evolution according to 

U(t) = D(a E (t))^*\ (B.5) 
where D{ole) is the displacement operator by an amount given by 

a E (t) = -Sp. f At'E(t')e wt ' (B.6) 
a Jo 

in the interaction picture. Starting from an initial state |</>), the resulting state is 

D(a E (t))e^ \<t>) . (B.7) 
The overlap with the initial state is 

e"W(<f>\D( aE (t))\<l>), (B.8) 
for which the modulus squared gives the probability of overlap 

Po(t) = \{ ( f ) \D(a E (t))\ ( f ) )\ 2 . (B.9) 
For a squeezed vacuum state 

E(0) = E(0 10} = e« at2 -^ a2)/2 |0) , (B.10) 
one can use [67] 

E(— ^)a^S(^) = a^coshr — ae _z<|) sinhr 

E(-0aE(0 = acoshr+ tie 1 * sink r (B.ll) 

to find that 

Po(t) = exp {-\a E f^^ - (a|e» + (a^f^l) . (B.12) 
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The loss of overlap with the initial state will be the ensemble average of this 
quantity. In order to evaluate it, we need to find 

V^WI 2 ), (B.13) 



with the exponent given by 

\a E {t)\ 2 = ^4- f dt' [ t dt ,, E(t , )E{t ,, )e- iuj ^ t,, - t, \ (B.14) 

h Jo Jo 

We can relate the expectation value of the exponential in equation B.13 to the 
exponential of the expectation value of its argument if we assume that the averaging 
time is short compared to the time scale over which the level populations vary, but 
large compared to the correlation time of the fluctuations. The expectation value of the 
argument can be found to be [68] 

<* E (t)\ 2 ) = ^S E {u z )t, (B.15) 

where Se(u z ) is the spectral noise density in the trap (see equation [I]) . 
This gives a decay constant (considering only the \cie\ 2 term): 
, g 2 cosh(2r) cosh(2r) 

T = 4^ ?b(w * )_ 2~ = IW_ 2~ ' (R16) 
with Tq^i defined by equation [TJ 



The exponent in equation B.12 also contains two terms with o? E and (cV^) 2 . The 
a 2 E term contains the integral 

f dt' I' dt"E(t f )E(t")e- iuJz V' +t '\ (B.17) 
Jo Jo 

This is a symmetric integral, in which we can make the substitution r = t" — t f to obtain 

* dtV w drE(t f )E(t r + r)e- ZUJzT . (B.18) 
J-t' 



Note the complex exponential in the integral over not present in equation |B.14[ Under 
the same approximation as before, and again taking the ensemble average, we obtain 
an additional sine function, meaning that this term tends to zero for t 3> 2tt/uj. This 
reflects the fact that the phase of the field in the Hamiltonian is uncertain and will 
be distributed over the full range from zero to 2tt. Hence, it will be averaged to zero 
over many realizations or over long interaction times. Therefore, taking the ensemble 
average, we see that the loss of probability of overlap is at a rate r _1 as given in equation 
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